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ABSTRACT 

Measurements  of  Pu-Be  neutron  age  to  indium  resonance 
and  thermal  migration  area  were  made  in  pure  water  at  2k °C 
and  in  mixtures  containing  1.16,  8.96-,  27.81,  ^6.65, 
70.09,  and  80.09  grams/liter  of  boric  acid  at  40°C.   Thermal 
diffusion  lengths  were  calculated  from  these  measurements 
and  the  dependence  on  poison  concentration  observed.   It 
was  concluded,  both  from  failure  of  the  age  measurement 
and  the  irresolute  thermal  flux  behavior  at  the  heaviest 
poison  concentration,  that  the  "maximum  absorption"  limit 
was  exceeded  in  this  experiment.   The  absorption  cross- 
section  of  the  supporting  medium  corresponding  to  the  80.09 
g/1  concentration  was  0.6120  cm""  .   The  present  results, 
together  with  previous  experimental  data,  provide  an  esti- 
mate of  diffusion  length  at  the  absorption  limit  of  0.5  cm, 
slightly  less  than  that  corresponding  to  the  theoretical 
Corngold  limit  (0.6  cm). 
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I.   INTRODUCTION 

It  has  been  shown  analytically  by  Corngold  [1]  that  if 
the  absorption  concentration  of  a  supporting  medium  is 
increased  beyond  a  certain  value  (assuming  1/v  absorption 
in  an  infinite  medium),  no  thermal  neutron  diffusion  length 
exists.   The  practical  implication  of  this  is  that  under 
such  extreme  circumstances  the  spatial  distribution  of 
thermal  neutrons  cannot  be  resolved.   Table  1  lists  theo- 
retical values  of  the  minimum  diffusion  length  which 
corresponds  to  the  "maximum  absorption"  limit  for  the  more 
common  moderators  [2]  . 

TABLE  1 


Moderator  Minimum  Diffusion 

Length  (cm) 


C  23 

Be  18 


H20  0.6 

D20  2.8 


The  dependence  of  thermal  neutron  diffusion  length  upon 
the  absorption  concentration  in  water  has  been  the  subject 
of  several  experiments  [33  4,  51 .   In  all  of  these  experi- 
ments the  limiting  value  of  diffusion  length  is  assumed, 
but  not  proved. 


The  most  extensive  measurements  of  diffusion  length  in 
water  poisoned  with  a  1/v  absorber  have  been  those  by 
Goddard  and  Johnson  [3].   Their  work  included  measurements 
in  water  containing  dissolved  boron  which  was  ideally 
suited  for  the  purpose.   Boron  is  a  strong  1/v  absorber, 
i.e., its  cross-section  varies  inversely  with  neutron  velo- 
city, or  as  the  square  root  of  neutron  kinetic  energy. 
Therefore,  the  homogenious  mixing  of  small  amounts  with 
pure  water  has  negligible  effect  on  the  average  scattering 
cross-section  but  drastically  increases  the  absorption 
cross-section . 

Goddard  and  Johnson  experimented  with  a  wide  range  of 
boron  concentrations  at  temperatures  of  22.3°  and  65°C. 
Their  most  heavily  absorbing  solution  contained  69.3  g/1 
of  boric  acid  (H,BCO  and  gave  a  diffusion  length  measure- 
ment of  O.658I  cm  (within  the  Corngold  limit).   Similar 
measurements  using  dissolved  cadmium  at  concentrations  of 
55.5  and  82.2  g/1  resulted  in  measurements  of  0.5983  and 
0.4925  cm  respectively,  both  of  which  lay  outside  the 
Corngold  limit.   They  did  not  report  any  difficulties  in 
obtaining  these  latter  diffusion  lengths. 

This  research  project  aimed  at  experimentally  determin- 
ing the  "maximum  absorption"  limit  by  increasing  the  boron 
concentration  in  pure  water  beyond  that  attained  by  Goddard 
and  Johnson  until  the  diffusion  length  measurement  failed. 

This  experiment  was  conducted  using  a  plutonium-beryHium 
source  which  emitted  5 A   x  106  fast  neutrons  per  second. 


It  will  be  explained  In  the  following  sections  why  this 
fast  neutron  source  dictated  an  unconventional  and  indirect 
method  of  determining  diffusion  length.   In  the  process, 
information  concerning  the  dependence  of  neutron  migration 
area  and  age  on  the  poison  concentration  is  obtained  for 
the  first  time. 


II.   THEORY 

A.   TWO-GROUP  DIFFUSION  THEORY 

The  thermal  neutron  flux,  in  most  cases,  arises  from  a 
source  of  neutrons  having  energies  many  orders  of  magnitude 
larger  than  the  energy  of  neutrons  in  thermal  equilibrium 
with  the  medium  (about  .025  eV  at  25°C)  [6].   The  Pu-Ee 
(a,n)  source  used  in  this  work  yields  neutrons  with  an 
average  kinetic  energy  of  about  *J .  2  MeV  and  the  neutron 
spectrum  ranges  from  roughly  zero  to  about  10.6  MeV.   The 
calculated  spectrum  is  shown  in  Figure  1.  [7] 
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Figure  1.   Calculated  Neutron  Spectrum 
for  a  Pu-Be  Source. 


The  analytical  model  which  describes  the  neutron  popula- 
tion as  two  separate  groups,  one  fast  and  one  thermal,  gives 
solutions  which  approximate  the  behavior  of  neutrons  in  real 
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assemblies  and  for  many  practical  calculations,  is  suffi- 
ciently accurate  [8].   In  the  two-group  treatment,  neutrons 
of  thermal  energy  constitute  one  group  while  all  those  of 
higher  energy  belong  to  the  fast  group.   Neutrons  are 
considered  removed  from  the  fast  group  by  scattering 
collisions  which  reduce  their  energy  to  the  thermal  region. 
Since  absorption  cross-sections  of  moderators  for  high 
energy  neutrons  are  very  small,  any  loss  of  the  fast 
neutrons  by  absorption  is  neglected. 

The  disappearance  of  neutrons  from  the  fast  group  is 
described  in  terms  of  a  "slowing  down"  cross  section,  Ep 

16]  ,   By  denoting  the  fast  neutron  flux  by  <J>p  (neutrons/ 

2 
cm  -sec),  Ep  is  defined  such  that  E„4>„(r)  gives  the  number 

of  neutrons  per  cnrVsec.  which  slow  down  out  of  the  fast 

group  at  the  point  r  cm  from  the  source. 

If  an  isotropic  point  source  at  r  =  0  emits  S  fast 

neutrons/sec  in  an  infinite  moderator,  the  fast  flux  at 

r  f1   0  can  be  described  by  [7] 

DpV2<t>p(r)  -  Epcf>p(r)  =  0  (1) 

where  V  4>(r)  is  the  Laplacian  of  the  flux  and  Dp  is  the 
diffusion  coefficient  for  the  fast  neutrons.   The  source 
at  r  =  0  is  taken  into  account  by  means  of  a  source  con- 
dition which  will  be  stated  later  in  this  section. 

The  quantity  Dp  can  be  derived  by  considering  Fick's 
Law  of  Diffusion  which  states  that  the  neutron  current 
density  is  proportional  to  the  negative  gradient  of  the 
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flux.   The  proportionality  constant  has  dimensions  in 

length  and  is  called  the  diffusion  coefficient,  denoted 

by  D  in  the  more  general  treatment  of  neutron  diffusion. 

Accordingly,  the  fast  neutron  current  density  can  be 

represented  by  [6] . 

En 

Fast   Neutron    _    n  n,  #  %      r  u  rtrc.\r7/L/      x?\    xitt 

r>  +-  r^    •  *.      =  -  D-oVcp-clr )  =  -  /  D(E)V(p(r,E)  dE 

Current  Density      F  F        p  ' 

fcth 
where  the  integral  is  evaluated  over  the  energy  region  of 

the  fast  f roup ,  i.e.,  from  thermal,  E,,  ,  to  source  energy, 

Efi .   The  fast  flux  is  defined  by  the  integral 

E 
<fu(r)  =  /  °  4>(r,E)  dE 

and  therefore 

E0 

/   D(E)V(j)(r,E)dE 

D   =  -^ 

F     E 

/  °  V<Kr,E)  dE 


E 


th 


In  two-group  theory  it  is  assumed  that  the  variation  in  the 
neutron  flux  spectrum  from  point  to  point  in  the  medium  can 
be  neglected,  that  is  the  flux  is  completely  separable  in 
space  and  energy.   Thus, 
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E 
/  °  D(E)4»(E)dE 

dp  =  "F (2) 

Eo 

/    <|)(E)  dE 
th 

where  <j>(E)  is  the  energy  dependent  flux. 

Returning  to  Equation  (1),  and  using  the  boundary  con- 
dition that  4>p(r)  must  remain  finite  as  r-*°°,  the  general 
solution  can  be  obtained: 

-r//r^ 
♦p(r)  =  C  |  ^  (3) 

where  the  parameter  i™  is  known  as  the  neutron  age  to 
thermal  energies,  defined  here  as  [9] 

•  DP 
tt  ■  ~  (4) 

The  constant  C  is  found  from  the  source  condition  which  is 
obtained  by  drawing  a  small  sphere  around  the  source  and 
computing  the  number  of  neutrons  that  pass  through  its 


surface  per  second.   Considering  the  magnitude  of  the  cur- 

P 


2 
rent  density  in  fast  neutrons  per  cm  /sec,  J,-,(r),  and  a 


2 

sphere  of  radius  r,  this  number  is  ^Jur  Jp(r).   Then,  in 

the  limit,  as  r ->  0 ,  the  source  condition  is 

lim  W2  Jw(r)  =  S  (5) 

r+0       * 
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Thus 

-r//x 


d<f>    (r)         D„Ce  T 

v*>  - -* -fc- - -^ — '^=+1 

and 

S 


L/tt 


C  "  ^Fa 


Finally,  the  fast  flux  is  given  by 

Considering  now  the  thermal  neutron  group,  the  time- 
independent  diffusion  equation  which  describes  the  thermal 
flux  can  be  written 

DV2  c|>T(r)  -  Za  <J>T(r)  +  qT(r)  =  0  (7) 

where  cf>m(r)  is  the  thermal . flux,  D  is  the  average  value  of 

the  diffusion  coefficient  in  the  thermal  region,  Z   is 

a 

the  average  macroscopic  thermal  absorption  cross  section 
and  qm(r)  is  the  thermal  slowing  down  density.   The  quantity 
qT(r)  is  equivalent  to  the  production  rate  of  thermal 
neutrons  per  unit  volume  and,  for  the  case  of  a  point  source 
of  fast  neutrons  in  an  infinite  moderator,  is  equal  to 
E  cf)  (r),  the  number  of  neutrons  slowing  down  out  of  the  fast 
group . 

Equation(7) can  be  rewritten  as 

<t>m(r)    £„  4>  (r) 
V2  <Mr)  -  -2-p-  +  -^ =  0  (8) 

1  v  D 
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where  <J>-p(r)  is  given  by  Equation  (2)  and  LT 


.-]   is 


the  thermal  neutron  diffusion  length.   By  requiring  the 

neutron  flux  to  remain  finite,  the  homogeneous  solution 

becomes 

-r/L^ 
K,  e     -1 

Vr)H  "  -^ '  <9> 

Since  the  source  term  is  proportional  to  the  fast  flux,  a 
particular  solution  is  assumed  of  the  form 

-r//rT 
K?  e 
*T(r)p  =  ,  (10) 

By  substituting  this  expression  into  Equation  (8),  the 
constant  K?  is  found  to  be 

-S  L  2 
K   = -i .  (11) 

4ttD  (LJ  -  tt) 

The  general  solution  to  Equation  (8)  can  now  be  written 

-r/LT       2   -r//^ 
K,  e        S  Lm   e 

♦T(r)  =  -^ ±L p (12) 

1  r       ilTrr  D(LT^  -xT) 

where  the  constant  K-.  can  be  determined  by  a  source  con- 
dition as  in  the  fast  flux  case.   However,  since  the  source 
only  produces  fast  neutrons,  the  required  source  condition 
for  the  thermal  flux  is 

lim  llTrr2  JT(r)  =  0  (13) 

r+0 
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where  Jm(r)  is  the  magnitude  of  the  thermal  neutron  current 
density.   Thus, 

_  d<f>T(r) 

Jn,(r)  =  -D  i 

1  dr 

and  the    constant   K,    becomes 

S   L    2 

K      =    —^ ^ (14) 

IttD  (Lt^  -tt) 

Finally,    the   thermal    flux   is    given  by 

„  -r/L^        -v/yfTZ 

*rp(r)    =   3—% (e  l    -e  X    )         (15) 

Hirr   EaLT^(l-TT) 

Equation  (15)  describes  the  thermal  neutron  flux  in 

2 
neutrons  per  cm  /sec  in  an  infinite  moderator  using  two- 
group  theory . 

It  follows  from  the  development  of  Equation  (15)  that 
it  is  also  possible  to  describe  the  slowing  down  of  the 
fast  neutrons  and  their  subsequent  diffusion  as  thermal 
neutrons  as  two  distinct  processes.   The  neutron  age,  Trps 
is  related  to  the  distance  traveled  from  the  point  where 
a  fast  neutron  first  appears  to  the  point  where  it  slows 

down  to  thermal  energies.   Similarly,  the  square  of  the 

2 
thermal  neutron  diffusion  length,  LT  >    corresponds  to  the 

distance  traveled  from  the  point  where  a  neutron  becomes 

thermal  to  the  point  where  it  is  finally  absorbed.   A 

quantity  MT  ,  known  as  the  thermal  migration  area,  takes 
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into  account  the  total  neutron  travel  during  moderation 
and  diffusion.   Accordingly, 

MT2  =  LT2  +  tt  (16) 

where  MT  (migration  length),  /rT  (slowing  down  length)  and 
LT  all  have  dimensions  In  length. 

B.   ONE-GROUP  DIFFUSION  APPROXIMATION 

Although  the  two-group  theory  gives  more  accurate  flux 
solutions,  the  one-group  approximation  readily  accommodates 
the  "first  scattering  source"  theory  which  is  useful  for 
our  analysis.   In  one-group  diffusion  theory  all  neutrons 
are  lumped  together  into  one-group,  namely  thermal.   This 
approximation,  while  limiting  the  description  of  the  neutron 
population  to  spatial  dependence,  allows  attention  to  be 
focused  on  the  diffusion  process. 

In  the  two-group  model  of  the  previous  section,  neutrons 
entered  the  thermal  group  as  a  result  of  slowing  down  out 
of  the  fast  group,  and  the  thermal  slowing  down  density, 
qrp(r),  appeared  as  the  source  term  in  the  thermal  diffusion 
equation.   In  the  one-group  treatment,  however,  when  the 
actual  neutron  source  yields  fast  neutrons,  the  source  term 
used  for  the  diffusion  equation  must  be  altered  to  account 
for  the  behavior  of  the  fast  neutrons  prior  to  reaching 
thermal  energies.   In  the  particular  case  of  a  point  source 
emitting  S  fast  neutrons/sec  in  an  infinite  media,  the 
slowing  down  density  at  r  is  given  by  the  relation  [10] 
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q(r,T)  =  ^ ^7p  (17) 

where  t  is  the  neutron  age.   A  small  t  means  that  the  neu- 
trons have  suffered  little  slowing  down  and  so  they  have 
not  diffused  far  from  the  source.   Larger  values  of  x 
correspond  to  neutrons  which  have  been  slowing  down  longer 
and  therefore  have  had  an  opportunity  to  diffuse  further 
from  the  source.   When  the  neutron  age  to  thermal  energies, 
Tm,  is  substituted  into  Equation  (17),  the  resulting  expres- 
sion gives  the.  dispersion  of  thermal  neutrons  and  can  then 
be  used  as  the  source  term  for  the  diffusion  equation. 

Equation  (17)  is  based  on  the  Fermi  age  model  which 
assumes  that  the  slowing  down  process  is  a  continuous  one 
[9] .   In  hydrogeneous  media,  however,  the  slowing  down 
process  is  not  continuous  due  to  the  fact  that  a  neutron 
may  lose  all  of  its  energy  following  a  collision  with 
hydrogen.   In  addition,  the  hydrogen  scattering  cross 
section  decreases  rapidly  with  energy  above  about  10  keV. 
As  a  result,  the  fast  neutrons  found  far  from  the  source 
are  those  resulting  from  high  energy  neutrons  which  have 
not  suffered  a  collision.   When  these  neutrons  finally  do 
collide  they  will  be  moderated  without  moving  very  far  from 
the  point  of  collision.   Accordingly,  a  fair  representation 
of  the  source  of  thermal  neutrons  which  can  be  used  as  the 
source  term  in  the  diffusion  equation  is  the  "first 
scattering  source"  [10]  which  for  a  point  source  is  given  by 
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q(r)  =  2  Es  e   S   * 

This  equation  gives  the  density  of  first  scattering  colli- 
sions at  distance  r  from  an  isotropic  point  source  emitting 
S  fast  neutrons/sec  in  an  infinite  medium.   The  parameter 
T.      is  the  macroscopic  scattering  cross  section  of  the 

fast  source  neutrons.   The  reciprocal  of  £   is  known  as 

^  s 

the  relaxation  length,  b,  defined  as  the  distance  the 
neutron  density  is  reduced  by  a  factor  1/e  of  its  initial 
value  due  to  scattering.   Replacing  I     by  1/b ,  the  source 
term  can  be  written 

K   e-r/b 

qCr)  =  -^—2 (18) 

r 

g 

where  Kn  =  -n— ■  Z  ,  a  constant. 
1    4tt   s  ' 

At  large  r,  therefore,  the  principal  variation  in  the 
slowing  down  density  arises  from  a  simple  exponential 
decrease  combined  with  the  inverse  square  characteristic 
of  a  point  source.   Moreover,  q(r)  is  not  a  function  of 
age,  t    since  every  neutron  is  presumed  to  slow  down 
fully  at  its  first  collision. 

The  thermal  diffusion  equation  for  the  one-group  approx- 
imation can  now  be  written  as 

p        0  (r)    Kie*-r/b 
V2(J)rp(r)  -  -i-«-  +    _  p    =  0  (19) 

1  v         Dr 
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To  satisfy  the  boundary  condition  that  4>T(r)  go  to  zero  as 
r  -*■  °°,  the  homogeneous  solution  becomes 


<Mr)H 


-r/LT 
Ke (20) 


A  particular  solution  is  then  assumed  of  the  form 

v"P  =  V  +  J?-'  <21) 

By  substituting  this  expression  into  Equation  (19),  the 
constants  C1  and  C?   are  found  to  be 


c1  =  0 


and 


_  Klb2LT  (22) 

2    T  2  .2 

The  general  solution  to  Equation  (19)  can  now  be  written 

a    i    \        K  e        1    T (23) 

♦«("-~ Dr2(LT2-b2)- 

If  the  diffusion  length,  LT,  is  larger  than  the  relaxa- 
tion length,  b,  the  first  term  of  Equation  (23)  dominates 
at  sufficiently  large  r,  the  thermal  flux  decays  expo- 
nentially, and  LT  can  be  evaluated  in  a  simple  routine 
manner.   This  is  the  conventional  method  of  determining  LT 

The  relaxation  length,  however,  depends  on  the  initial 
energies  of  the  neutrons  generated  by  the  source.   Fast 
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neutrons  from  a  Pu-Be  source  have  relaxation  lengths  of 
about  10  cm.  in  water  [11]  which  exceed  2.78  cm. ,  the  cor- 
responding diffusion  length  of  thermal  neutrons  [6] .   Con- 
squently,  the  thermal  flux  decays  in  a  more  complicated  way 
and  the  simple  theory  for  the  case  of  a  completely  thermal- 
ized  neutron  field  does  not  apply.   The  determination  of  L^ 
in  water  can,  however,  still  be  accomplished  by  measuring 
both  the  migration  area  and  the  neutron  age.  From  Equation 
(16)  the  diffusion  length  is  then  obtained  from  the  relation 

LT  =  (MT2  -  TT)3s  (24) 

C,   DETERMINATION  OF  NEUTRON  AGE 

The  neutron  age,  x,  between  two  energies  is  proportional 
to  the  mean  square  distance  travelled  by  a  neutron  as  it  is 

being  slowed  down  from  the  higher  to  the  lower  energy. 

_  2 
Denoting  this  distance  by  r„  ,  the  second  spatial  moment  of 

the  slowing  down  density  can  be  written  as  [9] 

9  /  r2q(r,T(En,E))dv 

V  (VE)  =   /q(r,T(Eo:E))dv   =  6t(VE)      (25) 

where  En  is  the  energy  of  the  neutron  at  its  point  of  birth 
as  a  fast  neutron  and  E  is  a  lower  energy  attained  as  a 
result  of  the  slowing  down  process;  dv  is  the  volume  element, 
and  q  is  given  by  Equation  (17).   Neutron  age  is  now  defined 
by  the  relation  [93 

t(Eq,E)  =  -|-  (EQ,E)  (26) 

Equation  (25)  applies  to  neutrons  of  any  energy  and, 
hence,  of  any  age.   In  actual  measurements,  however,  the 
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terminal  age  is  fixed  by  the  resonance  energy  of  the  detect- 
ing foil.   Consequently,  in  determining  the  age  to  thermal 
energies,  u,  the  spatial  distribution  of  the  slowing  down 
density  at  some  energy  above  thermal  is  obtained  and  then  a 
correction  is  applied  to  account  for  the  additional  age 
required  to  reach  the  thermal  region. 

The  slowing  down  density  of  source  neutrons  at  1.45  eV 
can  be  determined  by  activating  cadmium-covered  indium  foils 

placed  at  various  distances  from  the  source.   The  principal 

115 
indium  isotope,  IN    ,  has  a  strong,  isolated  absorption 

resonance  at  1.45  eV  while  the  absorption  cross-section  of 

cadmium  is  large  below  1.45  eV,  but  small  at  higher  energies 

The  cadmium  cover,  therefore,  prevents  thermal  neutrons 

from  reaching  the  indium  foil;  whereas  the  large  indium 

resonance  at  1.45  eV  is  primarily  responsible  for  any 

neutron  induced  activation, 

The  saturation  activity  of  a  cadmium  covered  indium 

foil,  ApD,  is  equal  to  the  rate  of  production  of  IN 

nuclei  and  can  be  written  as  [9] 

AnT)   =  V  /  <Kr,E)  £Q   (E)  dE  (2?) 

CD  aIN 

where  V  is  the  volume  of  the  foil,  Z    is  the  macroscopic 

115    IN 
absorption  cross-section  of  IN    ,  and  the  integral  Is 

evaluated  across  the  1.45  eV  resonance.   The  number  of 

neutrons  per  cm  /sec  which  slow  down  past  a  given  energy, 

E,  in  the  diffusion  approximation,  is  given  by 

q(r,E)  =  <|>(r,EK  E(E)E  (28) 

where  <f)(r,E)  is  the  energy  dependent  neutron  flux,  £  (E) 

is  the  macroscopic  scattering  cross-section,  and   £  is  the 


?p 


average  logarithmic  energy  decrement.   From  Equation  (28) 
it  follows  that 


ACD  "  V  ' 


rq(r,E)  ZaiN(E)-] 
5  VE) 


TT  (29) 


Assuming  that  q(r,E),  £,  and  £  (E)  are  constant  over  the 
small  energy  range  of  the  resonance  peak,  this  can  be 
written 


*cd  ■  vq<y-*5)  / 1    (e)  f         •      oo) 

s   s         IN 


That  is,  the  slowing  down  density  at  1.45  eV  is  proportional 
to  the  experimentally  measured  activity  of  the  foil.   By 
substituting  this  result  into  Equation  (25),  and  noting 
that  dv  =  llTrr  dr,  the  expression  for  the  age  to  indium 
resonance  in  spherical  coordinates  is 

°°  4 

r  2      '  r  ACD  dr 

tin  ~  "t~  "  6  ~~ -   :  ui; 

;0   r  ACDdr 

Thus,  by  knowing  only  the  relative  values  of  q(r,  1.45), 
the  age  to  1.45  eV  can  be  obtained  by  numerically  evaluating 
the  integrals  in  Equation  (31).   In  this  analysis  the  slowing 
down  density  is  assumed  to  be  approximated  by  Equation  (18). 

Another  method  for  measuring  neutron  age  [6] ,  which 
provides  only  a  rough  estimate  in  hydrogeneous  media,  is 
based  on  the  spatial  distribution  of  neutrons  of  a  given 
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energy  described  by  Equation  (17).   Taking  the  logarithms 
of  this  equation  and  considering  indium  resonance  energy 
leads  to 

£n  q(r)  =  const  -  r2/^TIN  .  (32) 

Since  the  resonance  activation  is  proportional  to  q(r) ,  a 
plot  of  in   Ap^  measured  at  various  points  from  the  source 
can  be  made.   For  moderators  which  satisfy  the  continuous 
slowing  down  theory  of  the  Fermi  Age  model  (e.g.  graphite) 
the  experimental  points  will  lie  on  a  straight  line.   For 
water,  however,  the  points  will  lie  above  a  straight  line 
except  in  a  very  narrow  band  about  20  cm.  from  the  source 
[12] .   This  is  due  to  the  unusual  nature  of  the  hydrogen 
scattering  cross-section  and  the  fact  that  individual 
neutron  energy  losses  are  so  large,  as  previously  discussed 
In  this  narrow  band,  however,  an  approximate  straight  line 
can  be  drawn  whose  negative  slope  will  provide  an  estimate 

°f  TIN« 

D,   DETERMINATION  OF  THERMAL  MIGRATION  AREA 

The  thermal  migration  area,  defined  by  Equation  (16), 
correlates  with  the  total  mean  square  distance  a  neutron 
travels  from  its  point  of  birth  to  its  point  of  capture. 
This  is  given  by  [11] 


f       Ag  r   dr 

?t2  -  -2= —  (33) 

/   AQ  r   dr 
0   b 


2k 


where  A„  is  the  total  saturated  activity  measured  as  a 
function  of  distance  r  from  a  point  source.   Equation  (33) 

can  be  obtained  in  a  manner  analogous  to  the  derivation 

_  2 
of  Equation  (31)  for  r„  .   Here,  however,  it  is  the  total 

relative  neutron  flux  which  must  be  known  to  evaluate  the 

integrals.   At  energies  that  activate  a  bare  indium  foil, 

this  total  flux  will  be  proportional  to  the  saturation 

activity.   Migration  area  can  then  be  determined  from  the 

definition 


-  2 
ft. 

lT     6 


Mm2  =  \-  (3M 
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III.   DESCRIPTION  OF  THE  EXPERIMENT 

Measurements  were  made  in  pure  H?0  at  24°C  and  in 
mixtures  containing  1.16,  8.96,  27.81,  46.65,  70.09  and 
80.09  grams  of  H  BO^  per  liter  of  HO  at  40°C.   The  last 
two  mixtures  contain  the  highest  boron  concentrations 
attained  to  date. 

The  experiments  were  conducted  in  a  cylindrical  plexi- 
glass tank,  1.^5  ft.  in  diameter  and  4.27  ft.  high,  which 
held  about  50  gallons  of  water.   This  tank  was  situated 
atop  a  smaller,  steel  tank  2.95  ft.  in  diameter  and  1.25 
ft.  high,  which  held  about  65  gallons  of  water.   The  base 
tank  served  to  contain  the  Pu-Be  source  and  also  provide  a 
protective  water  shield.   When  positioned  for  measurements 
in  the  base  tank,  the  source  was  0.625  inches  below  the 
base  of  the  plexiglass  tank  at  centerline.   A  polyethylene 
jacket  0.125  in.  thick  made  the  source  holder  watertight. 

A  single  steel  rod  immersed  in  the  tank  at  the  side 
heated  the  poisoned  mixtures  through  a  rheostat  control. 
This  heating  rod,  plus  a  continuous  running  electric  stirrer, 
maintained  the  medium  at  a  uniform  temperature  of  104°F 

(±0.5°). 

The  relative  flux  was  measured  with  bare  indium  foils, 
1.0  in.  in  diameter  and  .010  in.  thick.   Foils  used  in  any 
given  measurement  had  a  weight  spread  of  less  than  2.65%. 
Indium  resonance  activity  was  measured  with  the  same  foils 
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placed  in  cadmium  covers  0.025  in.  thick.   A  plexiglass 
foil  holder  arrangement  provided  accurate  spacing  of  the 
foils  in  the  medium  on  a  line  perpendicular  to  the  source. 
The  distance  between  adjacent  foils  was  0.50  inches  for 
all  measurements.   Figure  2  shows  the  experimental  appar- 
atus and  set-up. 

Nine  foil  positions  were  used  for  each  measurement. 
The  closest  measuring  position  to  the  source  was  8.26  cm. 
while  that  furthest  away  was  18.^1  cm.   Counts  of  the  cad- 
mium covered  foils  became  insignificant  beyond  this  last 
position . 

The  cadmium-covered  and  bare  indium  foils  were  irradi- 
ated separately  during  each  experiment  for  a  minimum  of  8 
hours.   After  irradiation  the  indium  was  allowed  to  decay 
for  3  minutes  to  avoid  interference  with  the  desired 
5^-minute  activity  from  the  13  second  half  life  of  IN 
[11] .   The  activity  of  the  foils  was  then  counted  from  each 
side  for  10  minutes;  the  relatively  long  counting  time 
having  been  chosen  to  improve  statistics.   The  counting 
system  consisted  of  an  end-window  Geiger  Mueller  tube  con- 
tained in  a  lead  shielded  housing,  and  an  electronic 
scaling  unit.   This  apparatus  is  summarized  in  Table  2 
and  pictured  in  Figure  3- 

Counter  background  was  measured  repeatedly  during  each 
run  as  an  erratic  line  voltage  perturbation  appeared  at 
infrequent  intervals.   Only  rarely  was  the  background 
stable  within  statistics  for  an  entire  run.   The  resulting 
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Experimental  Apparatus 
Figure  2 
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TABLE  2 
DETECTION  EQUIPMENT 


Counter 

Activity  Detected 

Type 

Shape 

Counting  Gas 

Manufacturer 

Amplifier 
Type 

Sensitivity 
Manufacturer 
Model 

Sealer-Timer 

Manufacturer 
Model 


Beta  and  Gamma 

Organically-quenched  GM 

Cylindrical 

Helium 

Atomic  Accessories,  Inc. 

Non-overloading 
10  mV 

Hewlett  Packard 
555^A 

Hewlett  Packard 
5201L 


necessity  for  interpolating  between  background  determina- 
tions represents  one  important  possibility  for  nonrandom 
errors  in  the  results. 

Prior  to  each  poison  experiment  samples  of  the  solution 
were  chemically  analyzed  for  boric  acid  content.   The  method 
of  analysis  was  boric  acid  titration  with  glycerine  which 
is  outlined  in  Appendix  A. 
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IV.   ANALYSIS  OF  DATA 

A.   DETERMINATION  OP  SATURATED  ACTIVITY 

Three  corrections  were  applied  to  the  foil  data.   The 
first  was  an  internormalization  based  on  foil  weight,  the 
second  corrected  for  absorption  of  indium  resonance  neutrons 
by  the  cadmium  cover;  and  the  third  corrected  for  counter 
resolving  time. 

The  weights  of  the  indium  foils  used  in  all  measurements 
varied  from  0.9282  to  0.9621  grams.   All  foil  weights  were 
therefore  normalized  to  a  weight  of  1.0  gram  and  the  result- 
ing corrections  applied  to  the  foil  counts.   The  highest 
value  of  the  weight  correction  factors  obtained  in  this 
manner  was  1.082. 

The  indium  resonance  counting-rate  correction  was 
determined  to  be  1.11.   This  value  was  obtained  from  a 
curve  of  correction  factor  vs.  indium  foil  thickness  con- 
structed for  a  cadmium  thickness  of  0.025  in.  as  outlined 
by  Tittle  [lH]  . 

A  number  of  determinations  of  counter  resolving  times 

-1  -1  c 

were  made  by  following  the  decay  of  5^  minute   IN    over 
three  hours.   An  average  value  of  *J00usec.  was  obtained  by 
solving 


A]_       A2 


1-Njt    1-N2t 


31 


where  A  is  the  extrapolated  activity,  N  is  the  mean  counting 
rate,  and  t  is  the  resolving  time.   A  correction  factor  of 
1/(1-Nt)  was  then  applied  to  all  counting  runs. 

The  effect  of  foil  activation  by  source  neutrons  above 
resonance  was  considered  negligible  as  the  closest  measure- 
ment point  was  greater  than  8  cm.  from  the  source  [12] . 
Also,  to  insure  that  the  close  spacing  between  foils  did 
not  contribute  error  due  to  shadowing,  activations  were 
made  in  pure  water  using  various  spacing  intervals. 

After  correcting  the  foil  count  as  described  above, 
the  saturated  activity  in  counts  per  minute  was  calculated 
from  the  equation 

_  /C 

S        -At,    -A(t0-t, )    -A(t0-t, ) 
•  (1-e    1)(e    2  1     -e     3  1    ) 

where  C  is  the  total  count  less  background,  A  is  the  disin- 
tegration constant  of  IN    ,  t,  is  the  irradiation  time, 
t?  -  t,  is  the  waiting  time  and  t^-t,  is  the  waiting  time 
plus  counting  time  1 13].   An  average  quantity  was  then 
obtained  from  the  values  of  A~  computed  for  each  foil  side. 
The  computer  program  used  in  this  calculation  is  included 
in  Appendix  F.   Appendix  I  provides  tables  of  the  computed 
data  for  all  runs. 

B.   DETERMINATION  OF  RELAXATION  LENGTH 

It  was  pointed  out  in  Section  II-B  that  in  hydrogeneous 
media,  the  slowing  down  density  at  large  distances  from  the 
source  can  be  approximated  by  the  expression 
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K   e"r/b 

q(r)  =  -J— ^ (18) 

r 


where  b  is  the  relaxation  length  of  the  fast  source  neutrons 
Since  in  the  framework  of  age-diffusion  theory,  the  flux 
is  proportional  to  the  slowing  down  density  at  a  fixed 
energy,  it  follows  that  the  cadmium-covered  activity  at 
large  distances  from  the  source  would  have  the  same  space 
dependence  described  by  Equation  (18),  namely 


K  P"r/b 


To  evaluate  b  from  the  corrected  data,  the  product 

2 

Ap^r  was  plotted  against  r  on  semi-log  paper.   This  plot 

showed  the  spherical  space  distribution  beginning  to 
decrease  exponentially  at  about  r  =  12.5  cm.   The  distribu- 
tion for  r  >  12.5  cm  was  then  fitted,  by  least  squares, 
with  functions  of  the  form  of  Equation  (35);  specifically 

Jin  ACDr  =  const  -  r/b  (36) 

The  distribution  for  r  >  12.5  cm  consisted  of  only  five  data 
points  which,  for  most  runs,  exhibited  large  statistical 
variation  due  primarily  to  the  low  foil  activity  in  this 
range.   Accordingly,  the  solution  resulting  from  fitting 
(36)  to  this  data  provided  only  a  rough  value  of  the  param- 
eter b.   A  more  accurate  determination  was  made  by  comparing 
ArD  distributions  calculated  from  Equation  (35)  (using 
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neighboring  values  of  b)  with  the  measured  distribution  on 

a  large  scale  semi-log  plot  of  A^  vs.  r.   Unlike  the 

2 
Ap^r   spherical  space  distribution  all  data  on  the  Arn  vs. 

r  plot  decreased  in  value  with  increasing  r;  this  permitted 
use  of  the  more  reliable  data  close  to  the  source  to  indi- 
cate  which,  if  any,  of  the  data  at  r  >  12.5  cm  were  incon- 
sistent with  the  decreasing  trend  of  the  overall  distribu- 
tion.  After  discarding  any  poor  data  indicated  by  this 
procedure,  the  calculated  A„n  distributions  were  normalized 
and  fitted  to  the  measured  distribution  for  r  >  12.5  cm 
until  the  closest  possible  fit  was  attained.   The  value  of 
b  used  in  calculating  the  distribution  found  to  give  the 

closest  fit  became,  then,  the  negative  reciprocal  of  the 

p 
slope  of  the  exponential  portion  of  theA^^r  vs.  r  curve. 

The  above  procedure  was  also  used  to  determine  relax- 
ation length  for  the  bare  foil  activity  (A„)  and,  thus,  the 

2 
slope  of  the  A„r   vs.  r  curve.   Foster,  in  measuring  the 

age  of  Na-Be  neutrons  [15] ,  assumed  no  loss  in  precision 

from  using  the  same  value  of  b  for  both  the  ACD  and  A„ 

distributions.   The  present  procedure,  however,  clearly 

showed  a  small  difference  in  the  two,  the  b  for  the 

resonance  distribution  being  less  in  every  case,  as  it 

should,   The  work  of  Rush  with  Ra-Be  neutrons  [16]  showed 

this  same  tendency.   Figure  h   is  a  plot  of  ACD  vs.  r  and 

A~  vs .  r  for  the  pure  water  run  which  includes  the  data 

computed  from  Equation  (35) 3    normalized  to  the  experimental 

data.   Figure  k   also  shows  the  marked  deviation  between 
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Figure    4 
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the  computed  and  experimental  data  for  points  near  the 
source.   This  is  to  be  expected  since  Equation  (35)  applies 
only  at  larger  distances  from  the  source.   Similar  plots 
for  the  poisoned  water  runs  are  presented  in  Appendix  B. 
The  values  of  b  determined  for  all  experimental  runs  are 
given  in  Table  3. 

TABLE  3 
RELAXATION  LENGTH 


Boric  acid 
(g/Jl) 

b( 
Ind 

A 

±0.03  cm) 
ium  Resonance 
ctlvation 

b( 
Tota 

±0.02  cm) 

1  Activation 

_  _  _ 

10.30 

10.50 

1.16  ±  .17 

10.25 

10.45 

8.96  ±.36 

10.25 

10.38 

27.81  ±  .79 

10.22 

10.36 

46.65  ±1.2 

10.19 

10.35 

70.09  ±1.3 

10.15 

10.35 

80.09  ±1.9 

10.15 

10.35 

C.   DETERMINATION  OF  AGE  AND  MIGRATION  AREA 

A  smoothed  curve  was  drawn  through  data  points  less  than 

2        2 
r  =  12.5  cm.  on  the  ACDr   and  Asr   versus  r  plots.   For 

r  >   12.5  cm.  the  data  were  extrapolated  to  infinity  by  means 

of  a  straight  line  having  a  slope  equal  to  the  negative 

reciprocal  of  the  parameter  b.   The  areas  under  the  smoothed 
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curves  were  found  by  integrating  the  straight-line  portions 
analytically  and  by  using  Simpson's  Rule  for  the  remainder. 
Points  from  these  curves  were  used  to  obtain  the  curves 
A„Dr   and  A„r  versus  r.   These  latter  curves  were  found 
to  approach  straight  lines  beyond  r  =  *J5  cm.  and  areas  were 

calculated  accordingly.   After  the  area  under  each  of  the 

_  2 
four  curves  had  been  determined,  the  second  moments  r^   and 

_  2 

r.   were  obtained  from  Equations  (25)  and  (33).   The  neutron 

age  to  indium  resonance  and  total  migration  area  then 

followed  from  Equations  (31)  and  (3*0.   The  computer  program 

used  in  these  computations  is  included  in  Appendix  G. 

2      2       4         I] 
Figures  5  and  6  are  plots  of  A„Dr  ,  A~r  ,  A^r  ,  and  A^r 

2        2 
versus  r  for  pure  water.   Plots  of  A^r   and  A„r   versus 

r  for  the  poisoned  experiments  are  given  in  Appendix  C. 

2 
Table  4  gives  the  values  of  ttn  and  M^   determined  for 

all  experiments.   The  errors  cited  have  been  estimated  from 

2        2 

A„Dr   and  A  r   versus  r  plots  by  using  the  individual, 

unaveraged  foil  points  to  construct  maximum  and  minimum 

2 
values  for  tjn  and  M„    [11] . 

A  rough  check  on  the  experimentally  determined  values 

of  age  was  made  by  following  the  procedure  discussed  in 

Section  II-D  with  Equation  (32).   Values  of  activity  as  a 

function  of  r  had  been  previously  computed  from  Equation 

(35)  in  determining  the  closest  value  of  relaxation  length 

to  fit  the  experimental  points.   This  data  was  plotted  on 

semilog  paper  against  r  .   A  narrow  band  between  r  =  18 

and  r  =  22  cm.  presented  the  closest  straight  line 
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Figure   6 
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TABLE  4 
MIGRATION  AREA  AND  NEUTRON  AGE 


Boric  Acid 

MT2(cm2) 

TIN(cm2) 

IN(cm  ) 

(Prom  Eq.  (32)) 

_  _  _ 

61.76  ±  1.84 

52.19  ± 0.79 

51.75 

1.16  ±  .17 

59.41 ±  0.68 

52.06  ±  2.34 

51.62 

8.96  ±  .36 

55.63  ±  1.17 

52.02  ±  3.18 

51.62 

27.81  ±  .79 

54.21  ±  2.89 

51.96  ±  3.06 

51.55 

46.65  ±1.2 

53.47  ±  3.13 

51.29  ±  3.42 

51.47 

70.09  ±1.3 

53.37  ± 1.57 

51.26  ±  3.11 

51.37 

80.09  ±1.9 

53.13  ±  3.13 

51.20  ±  4.21 

51.37 

TIN 


approximation.   The  slope  of  this  line,  equal  to  -  ^ 
from  Equation  (32),  was  then  determined  by  the  method  of 
least  squares.   A  plot  of  ACD  vs.  r2  for  pure  water  with 
b  =  10.3  is  provided  in  Appendix  D.   Table  4  includes  values 
of  TIN  obtalned  ln  tnis  manner  for  all  experimental  runs. 
The  least  square  machine  program  used  in  this  determination 
of  age  is  included  in  Appendix  H. 

D.   DETERMINATION  OF  DIFFUSION  LENGTH 

Prior  to  determining  LT  from  Equation  (24),  it  was 
necessary  to  apply  a  correction  to  tin  to  account  for  the 
additional  age  required  to  reach  the  thermal  energy  of 
0.025  eV. 
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Estimates  of  this  correction  for  hydrogeneous  media  are 

2  2 

found  in  the  literature  ranging  form  0 . 3  cm   [17]  to  1.8  cm 

[11] ,  the  latter  value  being  that  used  by  Valenti  for 

measurements  with  a  Pu-Be  neutron  source.   In  this  investi- 

gation  it  was  found  that  adding  a  correction  of  1.8  cm  to 

the  value  of  xTN  determined  for  pure  water  yielded  a  value 

of  2.786  for  L^.   This  is  in  excellent  agreement  with 

published  results  obtained  by  the  more  conventional  method 

of  measuring  the  exponential  attenuation  of  thermal  neutrons. 

It  was  also  noted  that  2.786  cm  was  in  precise  agreement 

with  the  value  of  L^  obtained  from  the  empirical  equation 

[14] 


LT    =  2.6*1  +  0.0061  T  cm  (37) 

H20 

for  a  temperature,  T,  of  24°C.   Accordingly,  a  correction 
of  1.8  cm  for  the  age  from  indium  resonance  to  thermal  ener- 
gies was  applied  to  the  experimentally  determined  value  of 
ttn  for  all  runs  and  L,p  was  then  computed  from  Equation  (2M). 
Results  are  presented  in  Table  5-   Errors  have  been  estimated 

from  a  combination  of  those  previously  cited  in  Table  k    for 

2 
tin  and  MT  . 

E.   DETERMINATION  OF  THE  THERMAL  FLUX  DISTRIBUTION  ACCORDING 
TO  TWO-GROUP  THEORY 

A  final  check  on  the  results  obtained  for  pure  H„0  and 
the  poison  mixtures  was  made  by  computing  the  thermal  flux 
distribution  from  Equation  (15).   The  experimentally  deter- 
mined values  of  LT  and  tt  were  used  in  the  equation  and  the 
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TABLE  5 

MEASURED  DIFFUSION  LENGTH  AT  2  4°C  IN 
PURE  WATER  AND  40°C  IN  POISONED  VJATER 


Boric  Acid  Lm(cm) 

g/£  -  1 


_  _  _  2.7862  ±  .14 

1.16  1  .17  2.3559  ±  .15 

8.96  ±  .36  1.3307  ±  .1^ 

27.81  ±  .79  0.6664  ±  .19 

46.65  ±  1.2  0.6213  ±  .20 

70.09  ±  1.3  0.5594  ±  .18 

80.09  ±  1.9  0.3640  ±  .30 


total  macroscopic  absorption  cross-section  of  the  poisoned 
medium  was  computed  from 

I     =    Ep  +  £a 

a    ahydrogen     Boron 

where  I  was  taken  to  be 

ahydrogen 

I  =  Z  +  NH  BO  X  ^°a 

ahydrogen    aH20     3   3     hydrogen 

and 

E       =  nh  RO  x  aa 
aBoron    H3B°3    aBoron 

The  parameter  a   is  the  microscopic  absorption  cross-section 

a  3 

and  N   DO   is  the  number  of  molecules  of  H3B03  per  cm 

H3B°3 
determined  from 
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PH  BO,,  A0 
m       =  2 2 

3   3   Molecular  weight 

where  p  is  the  density  in  grams/cm  and  A~  is  known  as 

23 

Avogadro ' s  number,  which  is  6.024  x  10   molecules  per  gram 

mole . 

The  computed  values  of  flux  were  then  normalized  to  the 

thermal  saturated  activity  which  was  obtained  by  subtracting 

the  covered  foil  from  the  bare  foil  data.   The  results  were 

plotted  together  on  semilog  paper  and  compared.   Figure  7 

gives  this  comparison  plot  for  pure  water.   Similar  plots 

for  other  runs  are  contained  in  Appendix  E.   Table  6  lists 

the  values  of  £   determined  for  each  run  together  with 

a  & 

the  maximum  deviations  found  between  the  experimental  data 
and  flux  distributions  obtained  from  Equation  (15). 

TABLE  6 

COMPARISON  OP  EXPERIMENTAL  RESULTS  WITH  SOLUTIONS 
OF  THE  TWO-GROUP  THERMAL  FLUX  EQUATION 


Boric  acid 

E  (cm"  ) 

a 

MAX  deviation 

_  _  _ 

0.0222 

11.7 

1.16  ±  .17 

0.0307 

14.2 

8.96  ±  .36 

0.0882 

12.5 

27.81  ±  . 79 

0.2270 

7.2 

46.65  ±  1.2 

0.3657 

3.5 

70.09  ±  1.3 

0.5384 

2.3 

80.09  ±1.9 

0.6120 

-  -  - 
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V.   COMPARISON  OF  RESULTS  WITH  THEORY  AND 
PREVIOUS  EXPERIMENTAL  MEASUREMENTS 

A.   AGE  AND  THERMAL  MIGRATION  AREA 

Theoretical  investigation  of  the  age  of  fission  neutrons 
in  pure  H?0  has  been  extensively  pursued  and,  for  this  case 
alone,  there  exist  calculations  of  the  actual  slowing  down 
distribution  which  can  be  compared  in  detail  with  measured 
distributions.   Comparable  results  have  not  been  obtained 
for  non-fission  sources  due,  primarily,  to  uncertainties  in 
the  energy  spectrums.   There  is  no  account  in  the  literature 
of  an  age  calculation  for  a  Pu-Re  source.   Coveyou  and 
Sullivan  [17] ,  however,  have  calculated  the  age  to  indium 
resonance,  ttn,  for  monoenergetic  sources  in  H?0  using  a 
Monte  Carlo'  age  code  which  provides  a  rough  comparison  for 
this  work.   The  plotted  results  of  their  calculations  demon- 
strate a  rapid  rise  in  age  with  increasing  source  energy; 
and  for  an  initial  neutron  energy  corresponding  to  the 

average  energy  of  a  Pu-Be  source  (4.2  MeV) ,  these  results 

2       2 
give  a  value  for  xTN  of  52  cm   (±.5  cm  )  . 

2 

Although  clearly  defined,  thermal  migration  area,  M^  , 

is  not  amenable  to  calculation  because  of  theoretical 
difficulties  encountered  in  describing  the  scattering  of 
thermal  neutrons  [18] .   As  was  borne  out  by  this  experiment, 
however,  thermal  diffusion  length  is  small  in  highly 

absorbing  systems  so  that  the  predominant  contribution  to 

p 

M    comes  from  the  slowing  down  process.   This  allows  some 


45 


2 
comparison  of  measurements  of  M„   and  calculations  of  T™ 

by  treating  the  thermal  diffusion  as  an  additive  correction. 

2 


Th 


e  only  measurements  of  MT   and  tt„  from  a  Pu-Be  source 


which  have  been  published  to  date  were  those  by  Valenti 
and  Sullivan  for  tin  [19],  and  by  Valenti  for  MT  [11]. 
Their  results  are  compared  with  those  obtained  in  this  work 
in  Table  7.   A  factor  which  accounts  for  some  of  the  dif- 
ference in  these  two  measurements  is  a  source  effect  cor- 
rection used  by  Valenti  and  Sullivan  [19] .   This  involved 
a  correction  to  the  indium  resonance  activity  to  account 
for  the  finite  size  of  the  source  and  was  intended  to 
approximate  a  distribution  closer  to  that  expected  from 
a  point  source.   Due  to  the  relatively  low  counting  rates 
encountered  in  this  work,  it  is  estimated  that  a  source 
effect  correction  would  have  increased  the  activity  of  the 
foil  nearest  the  source  by  less  than  1%    [11] . 

TABLE  7 

COMPARISON  OF  RESULTS  WITH  PREVIOUS  EXPERIMENTS 
Plutonium-Beryllium  Neutrons  in  Pure  Water 


Valenti  & 
Present  Results    Sullivan  [19]     Valenti  [11] 

Mm2(cm2)      61.76  ±  1.84  61.07 


T 
IN 


TXM(cm2)      52.19±0.79        52.8  ±2.0 


k6 


2 

In  Figure  8,  measured  value:;  of  the  parameters  M^   and 

ttn  are  plotted  as  a  function  of  poison  concentration.   The 
differences  between  the  two  curves  is  seen  to  reflect  the 

1/v  behavior  of  the  boron  absorption  cross-section.   At  low 

2 

thermal  energies  absorption  is  heavy  and  affects  only  M™  ; 

at  energies  above  thermal,  absorption  falls  off  rapidly  and 

2 
has  a  very  small  affect  on  both  M^  and  tt„. 

B.   THERMAL  DIFFUSION  LENGTH 

The  close  agreement  between  the  theoretical  value  of  L„ 
in  pure  HpO  at  room  temperature  with  that  obtained  in  this 
experiment  was  pointed  out  in  Section  IV-D. 

Previously  reported  investigations  of  LT  dependence  on 
the  poison  concentration  of  H?0  were  conducted  using  a  ther- 
mal reactor  as  a  neutron  source  [3,  1]  and  a  Sb-Be  (25keV) 
source  [5] .   As  a  result  there  existed  only  a  small  com- 
ponent of  fast-flux  in  the  region  of  measurement  which 
could  easily  be  eliminated  from  the  experimental  data. 
This  permitted  relatively  straightforward  measurements  of 
the  exponential  attenuation  of  thermal  neutrons  as  discussed 
in  Section  II-C  and,  consequently,  gave  final  results  of 
high  accuracy. 

Figure  9  compares  the  results  obtained  by  Goddard  and 
Johnson  [3]  with  values  of  L™  found  in  this  experiment.   The 
two  sets  of  data  are  found  to  be  in  good  agreement  when  the 
25°C  temperature  difference  is  taken  into  account.   Since 
their  results  were  obtained  at  65°C,  their  values  of  LT 
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should  be  higher  according  to  Equation  (37).   It  can  be  seen 
from  the  curve  representing  the  present  measurements  that 
absorbing  solutions  with  concentrations  of  70.09   and  80.09 
g/1  lay  outside  the  Corngold  limit  (0.6  cm). 

C.   TWO-GROUP  THERMAL  FLUX 

The  thermal  flux  computed  from  Equation  (15)  is  compared 
with  measured  values  of  saturated  activity  for  each  experi- 
mental run  in  Figure  7  and  Appendix  E.   Due  to  the  compli- 
cated manner  in  which  Pu-Be  neutrons  slow  down,  these 
computations  provided  the  only  method  of  examining  initial 
activity  data  for  agreement  with  a  distribution  predicted 
by  theory.   It  is  evident  from  the  maximum  deviations  given 
in  Table  6  that  the  measured  and  calculated  distributions 
come  closer-  in  agreement  with  increasing  poison  concentration. 

The  purpose  of  the  experimental  set  up  was  to  approxi- 
mate as  closely  as  possible  the  "infinite  medium"  condition 
which  was  a  basic  assumption  in  the  theoretical  treatment 
of  both  age  and  diffusion  length.   That  there  was,  in  fact, 
no  leakage  in  the  axial  direction  is  deemed  a  valid  conclu- 
sion from  the  absence  of  detectible  bare  foil  activity 
beyond  32  cm.  from  the  source.   An  indication  of  the  extent 
of  lateral  leakage  was  obtained  by  measuring  total  activity 
radially  outward  from  the  source  in  pure  HpO.   The  extrapo- 
lated distance  beyond  the  physical  dimensions  at  which  the 
flux  can  be  assumed  to  fall  to  zero  was  found  to  be  1.29  cm. 
This  result  is  displayed  in  Figure  10  from  which  it  can  be 
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10  15  20 

Radial  Distance  from  Source  (cm) 


22.86 


Figure  10 
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seen  that  even  in  pure  water  the  "infinite  medium"  assump- 
tion was  not  seriously  violated.   The  minor  differences 
between  the  measured  and  calculated  flux  distributions  for 
concentrations  of  ^6.65  and  70.09  g/1  indicate  the  better- 
ment of  the  assumption  at  the  heavier  absorber  concentration 
This  is  most  apparent  from  a  study  of  the  plots  in  Appendix 
E. 


RP 


VI.   CONCLUSIONS 

The  dependence  of  diffusion  length  on  poison  concentra- 
tion was  found  to  agree  well  with  that  reported  by  Goddard 
and  Johnson  131.   However,  because  diffusion  length  as 
obtained  in  this  experiment  is  the  square  root  of  a  small 
difference  between  two  large  numbers,  the  results  given  in 
Table  5  include  sizable  estimates  of  error.   Errors  in  the 
measurements  from  which  diffusion  length  was  obtained  were 
due  to  limited  experimental  conditions,  the  most  serious  of 
which  was  the  weak  neutron  source,  resulting  in  a  low  activity 
and  a  small  number  of  data  points.   This  in  turn  increases 
the  statistical  error. 

Previous  measurements  of  neutron  age  and  migration  area 
from  a  Pu-Be  source  in  pure  water  [11,  191  agree  with  the 
present  results  within  the  experimental  errors.   This  report 
is  the  first  of  measurements  of  these  parameters  as  a 
function  of  poison  concentration.   Accordingly,  the  results 
given  in  Table  4  should  be  of  considerable  value  in  extend- 
ing our  understanding  of  the  slowing  down  process.   The 
errors  associated  with  each  measurement  have  been  conserv- 
atively estimated  and  are  meant  to  serve  as  bounds  on  the 
error,  rather  than  a  probable  error. 

It  is  evident  from  the  experimental  data  plotted  for  the 
heaviest  concentration  (80  g/1)  that  it  was  impossible  to 
make  a  reliable  age  measurement  in  this  case.   The  value 
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given  for  this  run  is  therefore  merely  a  logical  estimate 
based  heavily  on  the  previous  measurements.   It  follows  from 
this  that  the  numerical  value  listed  for  diffusion  length 
at  this  concentration  is  also  unreliable.   However,  since 
it  can  be  assumed  that  each  of  these  values  approximates 
the  result  one  would  expect  to  obtain  had  a  measurement  been 
possible,  they  were  used  in  Equation  (15)  to  calculate  a 
thermal  flux  distribution.   The  result  of  this  calculation 
is  plotted  in  Figure  11  together  with  the  experimental  data 
through  which  no  curve  could  be  drawn.   Figure  11  vividly 
shows  the  irresolute  behavior  of  the  thermal  distribution 
at  the  80  g/1  boron  concentration.   It  is  concluded  from 
this,  and  from  the  failure  of  the  age  measurement  at  this 
concentration,  that  the  "maximum  absorption"  limit  was 
exceeded  in  this  experiment.   The  cadmium  measurements  of 
Goddard  and  Johnson  together  with  the  present  rerults 
indicate  that  the  diffusion  length  corresponding  to  this 
limit  is  about  0.5  cm,  slightly  less  than  the  theoretical 
value.   The  absorption  cross-section  of  the  supporting 
medium  at  the  time  the  final  measurements  were  attempted 
was  0.6120  cm"  . 
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APPENDIX  A 

METHOD  OF  CHEMICAL  ANALYSIS  FOR  DETERMINING 
BORIC  ACID  CONCENTRATION 


A.  PREPARATION  OF  STANDARD 

A  lr-gram  sample  of  H^BO^  (Mol  wt  61.84),  finely  ground, 
was  dissolved  in  300  ml  of  pure  H?0.  To  this  solution  in  a 
beaker  were  added  10  drops  of  phenolphthalein  indicator  and 
25  ml  of  neutral  glycerol.  From  a  buret,  standard  NAOH  was 
added  dropwise  until  a  distinct  reddish-pink  color  appeared. 
This  marked  the  true  end  point  for  boric  acid  and  established 
the  molarity  (Z)  of  the  standard  in  moles/liter. 

B.  ANALYSIS  OF  MODERATOR  CONCENTRATION 

To  25  ml  of  moderator  solution  were  added  indicator  and 
NAOH  as  described  above.   The  standard  alkali  was  then 
added  dropwise  until  the  true  end  point  was  reached. 

C .  EXAMPLE 

1.33  ml  of  NAOH  produced  end  point  in  solution 
containing  0.162  moles  H_BO~. 

Z  =  'qq^o  =  12.15  moles/liter   standard  alkali 

0.3  ml  of  standard  alkali  produced  end  point  in 
25  ml  of  sample  drawn  from  moderator 

(25  ml) (molarity  of  sample)  =  (.3  ml)(Z) 

molarity  of  sample  =  .1^6  moles/liter 

(61, 8^) ( .146)  =  9.03  g/l   HJBO,  concentration  of 
moderator 
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APPENDIX  B 

PLOTS  OP  INDIUM  RESONANCE  AND  TOTAL  ACTIVITY  VERSUS 
DISTANCE  PROM  THE  SOURCE  FOR  POISON  RUNS 
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APPENDIX  C 


PLOTS  OF  SPHERICAL  SPACE  DISTRIBUTION, 
ACDr2  AND  AQr2,  FOR  POISON  RUNS 


Figure  18 
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APPENDIX  D 

PLOT  OF  CALCULATED  INDIUM  RESONANCE  ACTIVITY 
VERSUS  r2  FOR  RELAXATION  LENGTH  OF  10.3  cm 


400 


1000 


2r       \2 
r  (cm) 


2000 


Figure  2k 


£o 


APPENDIX  E 
PLOTS  OF  THERMAL  FLUX  VS  r  FOR  POISON  RUNS 


Figure  25 


70 


!   '"■"!        _J_J_     _ 

_    }'T 

t   1. 

1 

__4_ 

T      1                     _ 

- 

t 

zj . 

.    __l       ! 

~i 

.   J  1__ 

2pM 

Measurec    Tf.erma 

r^/'ctitv 

xty-  - 

.-  »p- 

:g ~ 

'  BPI 

H_BC-    Goncer 

t ration 

-S.qfi 

:-t+t-: 
pr/1 

\  . 

H     ^         <! 

1 

^V^     ' 

-i     - 

i 

i  i 

l 

i 

\\» 

: 

V 

t 

I  -4- 



Xv 

\X  1 

-T         -.  vv 

\\ 

-  1—4- 

IC3 

^           H       """     4 

H 

Q     1-               ■     -4-t- 

Xi             i                    1 

in3 

rv~ 

_LU    ' 

ssfr 

>v 

lZ=s= 

X 

' 

th 

-}   —;              i-~i—t~t~- 

J— i — \r$~~  -            ^^m  B^ 

--  I  j  I 

(cpn) 

EHT      ill 

i4f+nrf 

mm 

SI 

:-  Caicu 

Two 
\        The 

lated 
Group) 

rmal    I 

1lux 

:^ttH 

4  — j— i — 1 \—  - 

1 — i „l_^ — , — j — 

SEE* 

*i 

. — »— 

-4= 

.:             3 

-f- — | 

— 

_           -       _! 

Is" 

\ 

t_ 

z[z 

v  ft 

3 

.1 

> 

-  X 

> 

X 

iSv 

v  X 

tV  r 

X 

:     zr: 

ZX 

s 

I 

V          X 

IS      v 

— 

-  ^ 

ln2 

■ 

fmT 

:Et      ~ 

10 

\ 

1 

■ 

-       ,      ,       — »— — 4—  — i— 

10 


12  llJ 

r(cm) 


16 


18 


Figure    26 


71 


Figure    27 
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APPENDIX  F 
COMPUTER  PROGRAM  1 

A.  PURPOSE 

1.  To  correct  counting  data  and  compute  the  saturated 
activity  in  counts  per  minute  (A-,-,  and  A,-,)  . 

2.  To  compute  the  approximate  mean  counting  rate  and, 
from  this,  determine  the  standard  deviation  based 

on  a  normal  distribution. 

2        2 

3.  To  compute  the  products  Ap~r   and  A~r  . 


B.   PROGRAM  NOTATION 


D 

T 

POS 

C,  CC 

ASAT,  BSAT 

AMEAN 

W 

SIGMA 

F 


rll6 


Decay  constant  for  IN' 

Counter  resolving  time 

Distance  from  source 

Ten  minute  counts  (side  1,  side  2) 

Saturated  activity  (side  1,  side  2) 

Approximate  mean  counting  rate 

Correction  for  foil  weight 

Standard  deviation 

Cadmium  correction  factor 


75 


p  APPENDIX   P 

COMPUTER   PROGRAM    1 


IMPLICIT    REAL*8(A-H,D-Z) 
WRITE( 6,1) 

1  FGRMAT(T5, 'DISTANCE' ,T25 , 'SAT    ACTI  VI  TY  »  ,T^6  ,  •  ME  AN'  , 
1T87, 'SIGMA'//) 

0=1.280-2 

T=6.CD-6 

D050C    M  =  l,2 

DG5C    N  =  l,9 

READ(5,2)    POS,C ,X,Y,Z,CC,YY,ZZ 

2  FORMAT( 6F9.0) 
READ(5,£)     W 

6  FORMAT(  1F9.0) 
DIST=POS*0. 25^001 
C  =  W*C/(  l.CDO-C*T) 
CC=W«CC/( 1.0D0-CC*T) 

ASAT=D*C/( ( 1.0DC-DEXP(-D*X) ) * ( DEXP ( -D*Y )-PEX P ( -0* I ) ) ) 
BSAT=OCC/(  (  1.0D0-DEXP(-D*X))*(DEXP(-D*YY)-CEXP(-D*ZZ) 
1)  ) 
SAT=(ASAT+BSAT)/2.000 
AMEAN=(C+CC) /2.0D01 
SIGMA=AMEAN**0.5D0 
F=1.11DC 
FSA7=F*SAT 
WRITE! 6,4)    DISTtFSATtAMPAN,SIGWA 

4  FORMAT ( T6,F5.2,6X,D2  5.16,6X,D25.16,6X,C2  5.16) 
A=FSAT*DIST**2.0D0 

B=FSAT*DIST**4.0D0 

8  FORMA T( /,T12,F5.2,6X,D25.16,6X,D25.16) 
WRITE (6,8)    DIST,A,B 

50    CON  TIN  IE 
D0  6C    N  =  l,9 
R£AP( 5,5)    POS,C,X,Y,Z,CC,YY,ZZ 

5  FORMAT( 8F9.0) 
READ(5,7)     W 

7  FORMAT(  1F9.0) 
0IST=P0S*0.254D01 
C=W*C/(1.0D0-C*T1 
CC=W*CC/(  1.0D0-CC*T) 

ASAT=D*C/( ( 1.0DC-DEXP(-D*X) ) * ( DEXP ( -0*Y )-DEXP ( -P* Z ) ) ) 
BSAT=D*CC/( ( 1.0D0-0EXP(-D*X) ) * ( DEX P( -D*YY ) -DEX P ( -D* ZZ ) 
1)  ) 
SAT=(ASAT+BSAT> /2.000 
AMEAN=(C+CC)/2.0D01 
SIGMA=AMEAN**0. 500 
WRITE( 6,3)    01  ST,  SAT,  A  n*F  AN,  SIGMA 

3  FORMAT ( T6, F 5. 2, 6X, 0  2 5. 16,6 X, 02 5. 15, 6 X, 02  5. 16) 
A=SAT*DIST**2.0D0 

B=SAT*DIST**4.0D0 
WRITE (6,9)    01  ST, A, B 

9  FORMAT(/,T12,F5.2,6X,D2  5.16,6X,P25.16) 
60    CONTINUE 

5CC    CONTINUE 
STOP 
END 
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APPENDIX  G 


COMPUTER  PROGRAM  2 


A.   PURPOSE 


1 .  To  compute  the  area  under  curved  portions  of  the 

2      2^         H 
curves  Ap„r  ,  A„r  ,  A„nr   and  A~r   versus  r  using 

Simpson's  Rule. 

2.  To  compute  area  under  straight  line  portions  of 

these  curves  analytically. 

2 

3.  To  compute  total  areas  and,  from  this,  M^   and 


T 


xTN  using  Equations  (13)  and  (18). 


*J .   To  compute  x^  by  correcting  tt„  for  additional  age 

required  to  reach  thermal  energies. 
5.   To  compute  diffusion  length  by  solving 

LT  =  (MT2  -  iT)h 


B.   PROGRAM  NOTATION 

2      2       4 
y(I)        Ordinate  values  of  A~Dr  ,  A~r  ,  A„Dr   and 

H  Abscissa  interval  lengths  in  cm.  used  to 

obtain  y(I)  values. 

N  Number  of  abscissa  intervals,  H. 

S(J)        Simpson  area  approximations 

B(I)        Relaxation  Length  values 

2     2      ^ 
Z(I)        Ordinate  values  of  A~Dr  ,  A~r  ,  A~,Dr 

and  A~r   corresponding  to  abscissa  values 
of  r  at  which  distributions  become  exponen- 
tial. 
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Ad) 

AA(I) 

MIG(I) 

AGE(I) 

Y(I) 

Area  under  exponential  portion  of  curves. 
Total  area  under  curves.  AA(I) =  S(J)  +  A(I) 
Thermal  Migration  area  values. 
Values  of  neutron  age  to  Indium  Resonance. 
Diffusion  Lengths. 
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APPENDIX  G 
COMPUTER  PROGRAM  2 

COMPUTING    SATURATED    ACTIVITY    TIMES    DISTANCE    SQUARED    AREAS 
IMPLICIT    REAL*8(  A-H,0-7) 
RE4L*eVIG 

D I  MEN'S  ION    Y(RO)  ,S( 23)  ,R(28)  ,7(28) , A(29) ,AA(23) , 
1MIG( 26) ,AGE(28) 
DO 80    J  =  1,1 4 
READ(5,10)     ( Y( I  )  ,1=1  ,11) 

10  FORMAT(  1F10.0) 
F=1.25 

N=10 

SUM2=C.CDC 

D0  5C    L=2,N,2 
5C    SUM2=     SLM2+Y(L) 

SUM3=0.CDC 

M  =  N-1 

D060  K=2,M,2 
6C  SUM3=SLM3+Y( K) 

S( J)=H/0.301*( Y( 1)+0.4D1*SUM2+0.2D1*SUM3+Y(N+1) ) 

WRITE(6,7C)     S(J) 
7C    F0RMAT(T19,  'SI KPSON    APPRCXI MATT CN= • f D23 . 16/// ) 
8C    CONTINUE 
COMPUTING    ACTIVITY    TIMES    DISTANCE    TO    THE    FOURTH    AREAS 

D090    J=15,28 

RtAD(5,ll)     ( Y( I )  ,1=1 ,19) 

11  FORMATt 1F10.0) 
H=2.5DC 

N=18 

SUM2=0.0D0 

D051    L=2,N,2 
51    SUM2=SIM2+Y(L) 

SUM2=0.CDC 

M  =  N-1 

00  61    K=3,",2 
61    SUM3=SLM3+Y( K) 

S( J)=H/C.3D1*( Y( 1 )+0.4Dl*SUM2+0.2Dl*SUM3+Y(N+l) ) 

WRITE(6,71)  S(J) 
71  FORMAT (Ti 9, • SI NP SON' ,025.16///) 
90  CONTINUE 
COMPUTING  AREAS  UNDER  LINEAR  PCRTICN  OF  CURVES 

READ(5,1)  (B(  I  )  ,Z( I  )  ,1=1,23) 

1  FORMAT( 2F1C.0) 
DO  2  1  =  1,25 

A(  I )=Z( I )/B( I  ) 
WRITE(6,3)  R( I  )  ,A(  I  ) 
3  FOkMAT(2Ei4.7///) 

2  CONTINL- 
00100  1=1,28 

100  AA( I )=A( I )+S(I ) 

D0200  1=1,13,2 

MIG( I )=AA(l+l4)/(AA(I)*6.0) 

WRITE (6,201)  MIG(I) 
201  F0RMAKT19,  »MIGRATICN  ARE  A=  »  ,D2  5  .1  6  ///  ) 
2CC  CONTINUE 

DG300  1=2,1^,2 

AGE(  I ) =AA(  1  +  14) /( AA(I )*6.C) 

WRITEJ  6,301)  AGE (I  ) 
301  F0RMAT(T19,«AGE=» ,025.16///) 
30C  CONTINLF 

D040C  I=2,lAf2 

Y(I)  =  (WIGU-1)-(AG3(I)+1.3))**0.5 

WRITE (6,  A  CI)  Y(I) 
4C1  FORMA! ( T19,  •DIFFUSION  LENGTH= • , 025 . 16//) 
4CC  CONTINUE 

STOP 

CND 
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APPENDIX  H 


COMPUTER  PROGRAM  3 


A.   PURPOSE 


1,  To  determine  slope  of  curve  ink,,     vs .  r   in  the 
range  18 .  ^4  £  r  <_   22  cm.  by  the  method  of  least 
squares . 

2.  To  compute  ttn  from  the  relation 


slope  =  - 


IN 


B,   PROGRAM  NOTATION 


F(I) 
G(I) 

X(I) 
Z(J) 


Relaxation  length  values 

Values  of  neutron  age  and  thermal  migration 
area 

Distances  from  the  source 

Solutions  to  equation 


Z  = 


Ce 


-r/b 


Y(I) 
SIMQ 


Natural  logarithm  of  Z  values 
Y  =  in   Z  =  a-,  +  a~  x 

Subroutine  for  solving  system  of  equations 
CA  =  B: 


i=l 


x. 

l 


7 

2   x. 

i=l   • 


T      x  ' 
i=l  1 


A 

— 

al 

a2 

B 


i  =  l 

7 

I 

i=l 


yi 


x.  y . 
i^i 


Slope  of  straight  line  solution  equal 
to  the  negative  reciprocal  of  ^tin 
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APPENDIX  H 
COMPUTER  PROGRAM  3 

DIMENSION  X(  9)  ,  Y(9)  ,7(P)  ,XC(2)  ,C(2,2) , 
1B(2)  ,A( 2)  ,D(7)  ,F( 7) ,G<9) 
REAO<5,5)  <D(I  ),F(I)  ,1=1,6) 

5  FGRMAT(2F9.Q) 
R£AD(5,5fc)  (G(I),I=1,6) 

5  6  FORMAT (F9.0) 
D06  L=l,6 
X(  1)=1£.41 
X( 2)=16.7 
X( 3)=19.2 
X  (  4  )  =  1 9  .  5 
X( 5)=19.7 
X( 6)=2C.O 
X(7)=22.C 
D09  J=l,7 

ZU)  =  (D(L)*EXP(-X(  J)/F(L)))/X<  J)**2.0 
9  CONTINUE 
DO  11  1=1,7 

11  Y(  I  )=ALOG<7.(  I)  ) 
D012  1=1,7 

12  X(  I  )  =  X( I )**2.0 
SUMMING  POWERS  OF  X(  I  ) 

DO  20  K=l,2 

XC(K  )  =  C.O 

CO  15     1  =  1,7 
15    XC(K)=XC(K)+X( I )**K 
20    CONTINUE 
SUMMING    Y( I ) 

YC  =  C0 

DO  2 2    1=1,7 
22    YC=YC+Y( I  ) 

YX=0.0 

D025     1=1,7 
25    YX=YX+Y( I )*X( I ) 
FORMING    TWO    EQUATIONS    CA=B    WITH    ELEMENTS     IN    COLUMNS 

D040  1=1,2 

D035  J=l,2 

JUGl=I+J-2 

IF( JUG1  )33,33,31 
33  C( 1,  1)=7.C 

GOTO  35 
31  C(  I, J)=XC( JUG1) 
35  CONTINUE 
40  CONTINUE 

B( 1)=YC 

B( 2)=YX 

CALL  SIMQ(C,B,2,KS) 

WRITE(6,5C)  KS 
50  FORMAT(T5,  »ERROR='  ,13) 

DO  5 2  1  =  1,2 
52  A( I )=R( I) 

hRITE( 6,54)  A(l) ,A(2) 

54  FORMAT (T5,  'A( 1) =  • ,E14.7/T5 , ' A( 2 ) =' ,E14.7) 
AGE=0.25/A(2) 

WRITE(6,55)  AGE 

55  F0RMAT(Ti5,,AGE=',E14,7//) 

6  CONTINUE 
STOP 

END 
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APPENDIX  I 
RESULTS  OF  COMPUTER  PROGRAM  1 

TABLE  8 
PURE  WATER;  INDIUM  RESONANCE  ACTIVATION 


r(cm) 

ACD(CPM) 

Mean  Count 

Standard 
Deviation 

8.26 

298.0691 

233.2133 

15.2713 

9.52 

204.1689 

157.3879 

12.5454 

10.79 

155.0969 

118.2727 

10.8753 

12.07 

124.0281 

89.9505 

9.4842 

13.34 

78.4197 

19.3193 

4.3954 

14.61 

67.7517 

50.4245 

7.1010 

15.88 

50.9003 

39.3977 

6.2767 

17.14 

42.2404 

32.7216 

5.7203 

18.41 

34.2155 

25.6226 

5.0618 

TABLE  9 
PURE  WATER;  TOTAL  ACTIVATION 


r(cm)      AQ(CPM)      Mean  Count     Standard 
S  Deviation 


8.26 

4459.2573 

3663.2972 

60.5252 

9.52 

3549.0267 

2918.2418 

54.0207 

10.79 

2747.3121 

2259.8505 

47.5378 

12.07 

1930.5774 

1664.8222 

40.8022 

13.34 

1437.4071 

1181.5408 

34.3735 

14.61 

1090.0382 

900.1418 

30.0023 

15.88 

852.3489 

692.3284 

26.3121 

17.14 

642.0987 

529.1137 

23.0024 

18.41 

541.2537 

449.0388 

21.1905 
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TABLE  10 
H-BO-  CONCENTRATION  1.16  g/1;  INDIUM  RESONANCE  ACTIVATION 


r(cm)      Apn(CPM)     Mean  Count     Standard 
CD  Deviation 


8.26 

503.1124 

379.3378 

19.4766 

9.52 

322.7215 

240.9674 

15.5231 

10.79 

233.0711 

176.0013 

13.2665 

12.07 

233.6401 

175.2742 

13.2391 

13.3^ 

172.9660 

129.3376 

11.3726 

14.61 

116.5811 

90.3156 

9.5035 

15.88 

74.4533 

57.2474 

7.5662 

17.14 

94.6393 

72.5835 

8.5196 

18,41 

48.3223 

36.5259 

6.0436 

TABLE  11 
H3B03  CONCENTRATION  1.16  g/1;  TOTAL  ACTIVATION 


r(cm)      AS(CPM)      Mean  Count     g™^n 


8.26 

6261.3213 

5172.9621 

71.9233 

9.52 

4148.3711 

3446.4026 

58.7061 

10.79 

3113.6634 

2596.8206 

50.9590 

12.07 

2229.4831 

1854.0849 

43.0591 

13.34 

1681.5223 

1400.5575 

37.4240 

14.61 

1265.8901 

1080.4300 

32.8698 

15.88 

986.9352 

842.4719 

29.0254 

17.14 

768.5511 

639.9204 

25.2966 

18.41 

635.0632 

527.7004 

22.9717 

83 


TABLE  12 
H  BO   CONCENTRATION  8.96  g/1;  INDIUM  RESONANCE  ACTIVATION 

r(cm)      Arn(CPM)     Mean  Count     Standard 
CD  Deviation 


8.26 

409.9467 

314.5537 

17.7356 

9.52 

291.4^77 

213.2623 

14.6035 

10.79 

229.5099 

167.0131 

12.9233 

12.07 

167.2182 

123.1536 

11.0975 

13.34 

142.2078 

103.2886 

10.1631 

14.61 

105.7302 

81.3866 

9.0215 

15.88 

67.2199 

51.4156 

7.1705 

17.14 

83.3672 

62.9565 

7.9345 

18.41 

44.4031 

33.5245 

5.7900 

TABLE  13 
H  BO^  CONCENTRATION '8.96  g/1;  TOTAL  ACTIVATION 


r(cra)      AQ(CPM)      Mean  Count     Standard 
S  Deviation 


8.26 

2713.8679 

2326.7305 

48.2361 

9.52 

1704.0191 

1413.4715 

37.5961 

10,79 

1188.5005 

986.2288 

31.4042 

12.07 

871.8115 

725.4009 

26.9333 

13.34 

639.1190 

530.0796 

23.0235 

14.61 

487.743 

405.6754 

20.1414 

15.88 

373.6443 

310.4774 

17.6204 

17.14 

279.7307 

226.8509 

15.0616 

18.41 

269.3002 

220.0322 

14.8334 
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TABLE  14 
H  B03  CONCENTRATION  27.81  g/1; INDIUM  RESONANCE  ACTIVATION 

r(cm)      Apn(CPM)     Mean  Count     Standard 
CD  Deviation 


8.26 

285.2747 

218.7960 

14.7918 

9.52 

236.3981 

183.6949 

13.5534 

10.79 

179.0414 

138.1119 

11.7521 

12.07 

135.2377 

103.4736 

10.1722 

13.3^ 

146.9174 

107.5823 

10.3722 

14.61 

89.3376 

66.5402 

8.5721 

15.88 

69.33^3 

51.0123 

7.1423 

17,14 

93.9194 

73.3271 

8.5631 

18. ill 

42.0807 

32.5463 

5.7049 

TABLE  15 
H^BO^  CONCENTRATION  27.81  g/1;  TOTAL  ACTIVATION 


r(cm)      AQ(CPM)      Mean  Count     Standard 
S  Deviation 


8.26 

1313.3395 

1118.5688 

33.4450 

9.52 

855.6882 

729.7254 

27.0134 

10.79 

658.1651 

558.1721 

23.6256 

12.07 

503.8277 

427.9465 

20.6868 

13.34 

360.5809 

306.2933 

17.5012 

14.61 

399.0305 

332.2687 

18.2282 

15.88 

206.7903 

176.8829 

13.2997 

17.14 

I86.5858 

158.2384 

12.5793 

18.41 

129.9547 

110.3157 

10.5031 
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TABLE  16 
H-BO-  CONCENTRATION  46.65  g/1;  INDIUM  RESONANCE  ACTIVATION 

r(cm)      A„n(CPM)     Mean  Count     Standard 
CD  Deviation 


8.26 

383.1435 

297.7103 

17.2543 

9.52 

282.4397 

212.2461 

14.5761 

10.79 

207.4628 

153.5804 

12.3927 

12.07 

162.0674 

119.9549 

10.9524 

13. 3^ 

150.3377 

113.8245 

10.6688 

14.61 

90.8321 

67.9496 

8.2432 

15.88 

95.1983 

71.8575 

8.4768 

17.14 

54.4682 

40.6000 

6.3718 

18.41 

79.3777 

57.6782 

7.5946 

TABLE  17 
H-BO^  CONCENTRATION  46.65  g/1;  TOTAL  ACTIVATION 


r(cm)      AS(CPM)      Mean  Count     ^TtTon 


8.26 

1006.1537 

835.3732 

28.9028 

9.52 

652.8737 

542.5276 

23.2922 

10.79 

455.5653 

377.8619 

19.4387 

12.07 

364.3542 

304.2721 

17.4434 

13.34 

310.0154 

251.2099 

15.8496 

14.61 

229.8385 

186.2102 

13.6458 

15.88 

199.8872 

159.1900 

12.6172 

17.14 

110.6291 

91.8081 

9.5816 

18.41 

88.0501 

73.0562 

8.5473 
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TABLE  18 
HJBCU  CONCENTRATION  70.09  g/1;  INDIUM  RESONANCE  ACTIVATION 

r(cm)      Apn(CPM)     Mean  Count     Standard 
CD  Deviation 


8.26 

379.0630 

285.6253 

16.9004 

9.52 

282.1943 

214.0142 

14.6292 

10.79 

230.4439 

165.9322 

12.8814 

12.07 

160.5295 

120.9588 

10.9981 

13.34 

106.3356 

79.5638 

8.9198 

14.61 

95.0191 

69.8891 

8.3599 

15.88 

72.4740 

55.6762 

7.4617 

17.14 

51.1528 

38.2288 

6.1829 

18.41 

54.6934 

41.7976 

6.4651 

TABLE  19 
H^BO^  CONCENTRATION  70.09  g/1;  TOTAL  ACTIVATION 


r(cm)      AQ(CPM)      Mean  Count     Standard 
S  Deviation 


8.26 

824.4368 

683.2489 

26.1391 

9.52 

558.2886 

462.7575 

21.5117 

10.79 

426.2032 

354.2170 

18.8206 

12.07 

316.5702 

261.6614 

16.1759 

13.34 

228.9763 

189.6963 

13.7730 

14.61 

192.6427 

159.1687 

12.6162 

15.88 

137.7566 

114.1718 

10.6851 

17.14 

93.6648 

79.3185 

8.9061 

18.41 

72.2664 

60.5687 

7.7826 
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TABLE  20 
H3B03  CONCENTRATION  80.09  g/1  J  INDIUM  RESONANCE  ACTIVATION 


r(cm)      Apn(CPM)     Mean  Count     Standard 
CD  Deviation 


8.26 

414.8513 

307.2195 

17.5277 

9.52 

256.6541 

196.2806 

14.0100 

10.79 

239.3121 

177.5645 

13.3253 

12.07 

161.0497 

122.6576 

11.0751 

13.31 

113.0485 

87.4346 

9.3506 

14.61 

99.8647 

77.9423 

8.8285 

15.88 

90.3185 

68.9324 

8.3026 

17.1^ 

60.4457 

45.0622 

6.7128 

18.41 

35.4877 

27.8769 

5.2799 

TABLE  21 
H-BO-,  CONCENTRATION  80.09  g/1;  TOTAL  ACTIVATION 


r(em)      AS(CPM)      Mean  Count     ^iaWon 


8.26 

743.8847 

614.4305 

24.7877 

9.52 

494.7427 

419.9927 

20.4937 

10.79 

360.5303 

309.5135 

17.5929 

12.07 

302.0725 

261.4078 

16.1681 

13.34 

212.3189 

175.3501 

13.2419 

14.61 

149.1620 

128.0910 

11.3177 

15.88 

127.8776 

105.3823 

10.2656 

17.14 

116.5993 

99.1096 

9.9554 

18.41 

70.7709 

53.4546 

7.3113 
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